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Abstract 

We report preliminary results on the decay B~ D^K*~ , B" DcpK*~ and its charge 
conjugate using a data sample of 95.8 million BB pairs recorded at the T(45') resonance with the 
Belle detector at the KEKB asymmetric e+e" storage ring. We find the branching fraction for 
B~ — > D^K*^ to be ^? = {5.2±0.5{stat) ±0.6(sys)) x 10~^ and the partial-rate charge asymmetries 
for DcpK*~ to be A = -0. 02±0. 33{stat) ±0. 07 (sys) and ^ = 0.19±0.50{stat)±0M{sys) 

where the indices 1 and 2 represent the CP=-|-1 and CP=— 1 eigenstates of the D° — system, 
respectively. 

PACS numbers: 13.25.Hw, 14.40.Nd 



*on leave from Fermi National Accelerator Laboratory, Batavia, Illinois 60510 
^on leave from Nova Gorica Polytechnic, Nova Gorica 



3 



The extraction of 03 [ij, an angle in one of the unitarity triangles of the Kobayashi- 
Maskawa (KM) quark mixing matrix P], is a challenging measurement even with modern 
high luminosity B factories. Recent theoretical work on B meson dynamics has demonstrated 
the direct accessibility of 03 using the process B~ — > DK*~ If the is reconstructed 
as a CP eigenstate, the b ^ c and h ^ u processes interfere. This interference leads to 
direct CP violation as well as a characteristic pattern of branching fractions. However, the 
branching fractions for D meson decay modes to CP eigenstates are only of order 1 %. Since 
CP violation through interference is expected to be small, a large number of B decays is 
needed to extract 03. Assuming the absence of — D'^ mixing, the observables sensitive 
to CP violation that are used to extract the angle 03 are 

^ B{B- Di^2K*-) - B{B+ Di^2K*+) 

2r sin 5' sin 03 
1 + + 2r cos 5' cos 03 




6 for Di 
5 + vr for D2 ' 



where Di and D2 are CP-even and CP-odd eigenstates of the neutral D meson, r denotes 
a ratio of amplitudes, r = \A{B~ D^K*~)/A{B~ D^K*~)\, and 5 is the strong 
interaction phase difference. The ratio r corresponds to the magnitude of CP asymmetry 
and is suppressed to the level of ~ 0.1 due to the CKM factor (~ 0.4) and a color suppression 
factor (~ 0.25). 

This measurement of the branching fraction for B^ D^K*~ and CP asymmetries for 
B^ — i> DcpK*~ is based on a data sample of 95.8 million BB pairs, collected with the Belle 
detector at the KEKB asymmetric-energy e+e" (3.5 on 8 GeV) coUider @ operating at the 
T(4S') resonance. At KEKB, the T(45') is produced with a Lorentz boost of (3'-f = 0.425 
nearly along the electron beamline. 

The Belle detector is a large- solid- angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50- layer central drift chamber (CDC), an array of 
silica aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of- flight 
scintillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) 
crystals located inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. 
The iron fiux-return located outside the coil is instrumented to detect mesons and to 
identify muons (KLM). The detector is described in detail elsewhere j^. 

We reconstruct mesons in the following decay channels. For the flavor specific mode 
(denoted by D/), we use D° K~ti~^, i^~7r+7r° and K~7i~^7i~7i~^ PJ. For CP =+1 modes, 
we use Di — > K~K^ and 7r"7r+ while, for CP = —1 modes, we use D2 — > K^tt^, K^ip, and 
K^uj. We reconstruct K*~ candidates from K^n' combinations where the candidate is 
formed from two oppositely charged pions having a vertex displaced from the interaction 
point in the direction of the momentum. The candidates are selected in the mass 
window of 0.492 GeV/c^ < M(7r+7r") < 0.505 GeV/c^. The candidate tracks are then 
kinematically constrained to its nominal mass value. The K*~ is required to have a mass 
within ±75 MeV/c^ of its nominal mass value. For the 7r° from the K^-k^-k^ decays, 

we require the 7r° momentum in the T(4S') center-of-mass(c.m.) frame be greater than 
0.2 GeV/c and the energy of each photon from the tt" be greater than 30 MeV. 

Well constrained reconstructed tracks that are not identified as electrons or muons are 
used as charged hadrons. For each charged track, information from the ACC, TOF and 
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specific ionization measurements from tlie CDC are used to determine a K/tt likeliliood 
ratio P{K/7t) = Lk/{Lk + -^tt), wliere Lk and Lj^ are kaon and pion likelilioods. For kaons 
(pions) from the K^tt^ mode we used the particle identification requirement of 

P{K/tx) > 0.4 (< 0.7). For kaons from the 1)° K-K+ mode we require P{K/tt) > 0.7 
while for pions from 7r^7r+ mode we require P{K/'k) < 0.7. 

The uj mesons are reconstructed from tt^t^^ti^ combinations in the mass window 
0.732 GeV/c^ < M(7r+7r"7r°) < 0.82 GeV/c^ with the charged pion particle identification 
requirement P{K/7r) < 0.8. To reduce the contribution from the non-resonant background, 
a helicity angle requirement |cos 6'hei| > 0.4 is applied, where 6'hei is the angle between the 
normal to the u decay plane in the u rest frame and the uj momentum in the frame. 
To remove the contribution from K*~p~^, we require the Kg7i~ invariant mass to be 

greater than 75 MeV/c^ from the K*^ nominal mass. 

The (j) mesons are reconstructed from two oppositely charged kaons in the mass window 
of 1.008 GeV/c2 < M{K+K') < 1.032 GeV/c^ with P{K/n) > 0.2. We also apply the 
helicity angle cut |cos 6'hei| > 0.4, where 6*1161 is the angle between one of the daughters in 
the (j) rest frame and the momentum in the frame. The candidates are required to 
have masses within ±2.5(T of their nominal masses, where a is the measured mass resolution 
that ranges from 5 MeV/c^ to 18 MeV/c^ depending on the decay channel. A mass and 
(wherever possible) vertex constrained fit is then performed on the remaining candidates. 

We combine the and K*~ candidates to form B candidates. The signal is identified 
by two kinematic variables: the energy difference AE = + E^*- — -E'beam and the beam 
energy constrained mass Mbc = \/ E^eum ~ \Pd Here E^ is the energy of the 

candidate, E^*- is the energy of the K*~ and -Ebeam is the beam energy, and ph and Pk*- 
are the momenta of the and K*~ candidates, respectively, all calculated in the cm. 
frame. Event candidates are accepted if they have 5.2 GeV/c^ < Mbc < 5.3 GeV/c^ and 
|Aii^| < 0.2 GeV. In case of multiple candidates from a single event, we choose the best 
candidate on the basis of a determined from the differences between the measured and 
nominal values of Md and Mk*-. Since B~ — > D^K*~ is a pseudoscalar to pseudoscalar- 
vector decay, the K*~ is polarized. We define the K*~ helicity angle cos 6'hei as the angle 
between one of the K*~ decay products in the K*~ rest frame and the K*~ momentum 
in the B rest frame. The K*~ helicity angle follows a cos^^hei distribution. We require 

I cos ^hell > 0.4. 

To suppress the large combinatorial background from the two-jet-like e+e~ ^ qq {<i = u, 
d, s or c) continuum processes, variables that characterize the event topology are used. We 
use a Fisher discriminant F, constructed from six modified Fox Wolfram moments |j7j and 
cosOb, the cosine of the angle of the B flight direction with respect to the beam axis, in a 
single likelihood ratio variable (LR) that distinguishes signal from continuum background. 
We apply a different requirement for each sub-mode based on the expected signal yield and 
the backgrounds in the Mbc sideband data. To give an example of the performance of this 
selection, the LR > 0.4 requirement keeps 89.4 % of the B~ — > D^[K~ K~^]K*~ signal while 
removing 69.2 % of the continuum background. 

The signal yields are extracted by a fit to the AE distribution in the region 5.27 GeV/c^ < 
Mbc < 5.29 GeV/c^. The signal is parameterized as a Gaussian with parameters determined 
from MC simulation. The continuum background function is modeled as a first order poly- 
nomial function with parameters determined from the AE distribution for the events in 
the sideband region 5.2 GeV/c^ < Mbc < 5.26 GeV/c^. Backgrounds from other B decays 
such as B- D*°K*-, where D*^ D^u^ and D*^ D^, are modeled as a smoothed 
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FIG. 1: The distributions in the B' D^K*- signal region for (a) AE with 5.27 GeV/c^ < 
Mbc < 5.29 GeV/c^. The soHd-hne shows the fit, the dashed-hne is the signal, the dot-dashed line 
is coninuum background and the hatched histogram shows the contribution from other B decays. 
In (b) we show Mbc distribution with a ztSo" AE cut where the solid-line shows the fit, the dotted 
line shows the continuum background and the dashed line is signal. In (c) we show the results of 
fits to the AE' distribution in bins of Ks'^~ invariant mass (points with error bars). In (d) we 
show the results of fits to the AE distribution in bins of cos ^hei (points with error bars). In both 
(c) and (d) the hatched histogram is a Monte Carlo (MC) simulation of B~ D^K*~ . 

histogram obtained from the MC simulation. The AE and Mbc distributions are shown in 
Figs. 1(a) and (b), respectively. We observe a signal of 169.5±15.4 events with 15. 4o" sta- 
tistical significance. The fit results are given in Table I. Statistically significant signals are 
observed for the B~ —>■ DfK*~ decay in all three decay channels. The final branching 
fraction is calculated from a weighted average of the results from these three channels. 

We remove the Kg7i~ invariant mass requirement and fit the AE distribution in bins 
of K^7i~ invariant mass and plot the signal yield from the fit for each bin. The result is 
shown in Fig. 1(c) and is consistent with a pure DK*~ signal MC simulation. We see no 
indication for non-resonant D^Kg7T~ production. Similarly, we remove the K*~ helicity 
angle requirement and fit the AE distribution in bins of cos ^hei- We then plot the signal 
yield from the AE fit for each bin. The result is shown in Fig. 1(d) and is consistent with 
the expectation for a pseudoscalar to pseudoscalar-vector decay. 

The following sources of systematic error are found to be sizeable: the tracking effi- 
ciency (1% per track), 7r° efficiency (4.8%), Kg efficiency (4.5%), fitting of the BB back- 
ground (5-7%) and particle identification (6-12 %). Other backgrounds including rare de- 
cays that could contribute to the AE signal region, are estimated from the sideband 
data (1 %). The uncertainty in the AE signal shape parametrization is determined by 
varying the mean and width of the signal Gaussian parameters within their errors. The 
uncertainty from the slope of the background is determined by changing its value by its 
error. Both of the resulting changes are included in the systematic error from fitting. The 



6 



combined systematic errors from these sources are 9.5% for B D'^lK 'n'~^]K* , 11.8% for 
B- D^[K-7i+7i^]K*- and 15.5% for 5" ^ D'^[K-7i+7i+7r-]K*- . The systematic error 
in the final branching fraction is calculated by weighting according to the efficiency times 
branching fraction of the three decay channels j8|. 




FIG. 2: The distributions in the DiK*^ signal region for (a) AE with 5.27 GeV/c^ < 

Mbc < 5.29 GeV/c^ (b) Mbc with a ±3a AE cut. The distributions in the D2K*~ signal 

region for (c) AE with 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ (d) Mbc with a ±3a AE cut. The 
solid-line shows the fit and the dashed-line shows the signal. 



TABLE I: The fit results for B DjK*^ decay. The signal yields, statistical significance, 
and branching fraction for each mode are given. The errors shown in the BF are statistical and 
systematic, respectively. 



decay mode 


Yield(A£; fit) 


Yield(Mfec fit) 


sig 


BR(10-^) 


Df K--K+ 


67.6±9.5 


75.1±9.5 


10.8 


5.0 ± 0.7 ± 0.5 


Df K-TT+TT^ 


66.3±9.6 


65.3±9.4 


9.5 


6.3 ± 0.9 ± 0.7 


Df — > K~Tr^7r^iT~ 


40.3±7.9 


43.0±8.1 


6.9 


4.3± 0.8 ± 0.7 


weighted mean 








5.2 ±0.5(stat) ± 0.6(sys) 



For B DcpK* decay the signal yields are extracted by a fit to the AE distribution 
in the region 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ where -0.2 GeV < AE < -0.1 GeV 
is excluded. The signal is parameterized as a Gaussian with parameters determined from 
MC simulation. The continuum background function is modeled as a first order polynomial 
function with parameters determined from the AE distribution for the events in the sideband 
region 5.2 GeV/c^ < Mbc < 5.26 GeV/c^. The AE and Mbc distributions are shown in 
Fig. 2. We observe a signal of 13.1±4.3 events for B~ — > DiK*~ and 7.2±3.6 events for 
B^ —>■ 1^2-^*^ with 4.3(7 and 2Aa statistical significances, respectively. The partial-rate 
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TABLE II: Yields, partial-rate charge asymmetries and 90 % C.L intervals for asymmetries. 



Mode 




N{B+) 


N{B-) 


Acv 


90 % C.L 


B^^ 


DfK*^ 


68.9±10.1 


95.3±11.3 


0.16 ± 0.09±0.08 


-0.04< Af <0.36 


B^^ 


DiK*^ 


6.7±3.0 


6.5±3.1 


-0.02 ± 0.33 ±0.07 


-0.57< Ai <0.53 


B^^ 


D2K*^ 


2.9±2.4 


4.3±2.7 


0.19 ± 0.50±0.04 


-0.63< A2 <1.00 



asymmetries Ai^2 are evaluated using signal yields obtained from separate fits to the 
and B~ samples. The results are given in Table II. The systematic uncertainty is from the 
intrinsic detector charge asymmetry (1 %), the B~ and B'^ yield extractions (4-7 %), and 
the asymmetry in particle identification efficiency of pions (1 %). The systematic error from 
yield extraction is calculated by changing the fitting parameters by ilo". At 90 % C.L 
intervals, we find 

-0.57 < ^1 < 0.53, 
-0.63 <A2< 1.00 

In summary, using 88 fb~^ of data collected with the Belle detector, we report a mea- 
surement of the exclusive decay mode B^ DfK*~ . This mode has been observed 
previously In this paper, we report a new measurement of the branching fraction 
B{B~ — >■ DfK*~). We also report the partial-rate charge asymmetries for the decay 
B^ DcpK*", where Dcp are the neutral D meson CP eigenstates. The measured 
partial-rate charge asymmetries Ai^2 are consistent with zero. 
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